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ABSTRACT 

We present projected velocity dispersion profiles for the Draco and Ursa Minor (UMi) dwarf 
spheroidal galaxies based on 207 and 162 discrete stellar velocities, respectively. Both profiles show a 
sharp decline in the velocity dispersion outside ^ 30' (Draco) and ~ 40' (UMi). New, deep photometry 
of Draco reveals a break in the light profile at ~ 25'. These data imply the existence of a kinematically 
cold population in the outer parts of both galaxies. Possible explanations of both the photometric 
and kinematic data in terms of both equilibrium and non-equilibrium models are discussed in detail. 
We conclude that these data challenge the picture of dSphs as simple, isolated stellar systems. 
Subject headings: dark matter — galaxies: individual (Draco dSph, Ursa Minor dSph) — galaxies: kine- 
matics and dynamics — Local Group — stellar dynamics 



1. INTRODUCTION 

The dark matter dominated Local Group dwarf 
spheroidals (dSphs) have emerged as valuable labora- 
tories in which to test dark matter models. Recently, 
the projected velocity dispersion profiles of the For- 
nax and Drac o dSphs have been obtained l)Matedll997t 
iKlevna et al.l |2001). Detailed modelling of the dis- 
crete stellar velocities in Draco enabled the hypothe- 
sis that mass follows light to be discarded at the 2. So- 
level and suggested that the halo density p{r) falls off 
more slowly than the light distribution with p(r) ~ 
r~^-^ (|Klevna et aL 2002). In this Letter, we present new 
observations of the Draco and Ursa Minor (UMi) dSphs 
which yield the velocity dispersion profiles of both galax- 
ies to the edge of their light distributions. The new data 
suggest the existence of kinematically cold populations in 
the outer parts of both galaxies. These data make it pos- 
sible to test the validity of isolated, equilibrium models 
of dSphs. 

2. OBSERVATIONS 

2.1. Discrete Velocities 

We observed the Draco and UMi dSphs with the mul- 
tifibre instrument AF2/WYFF0S on the William Her- 
schel Telescope on La Palma on 20-23 June 2003 and 6-11 
May 2003, respectively. We drew our Draco targets from 
the Sloan Digital Sky Survey, while we took our UMi tar- 
gets from our own KPNO 4m MOSAIC imaging. In each 
case, we identified potential targets by drawing a polygon 
around the giant branch of a V, ^ — / color- magnitude di- 
agram with a faint magnitude limit of T^ < 20. The data 
were reduced using the WYFFOS-specific WYFRED data 
reduction package in TRAP and cross-correlated with the 
two blue-most lines of t he Calcium tr i plet, in the same 
manner as described bv IKlevna et alJ l)2002l ^. To deter- 
mine the final member list for each dSph we assumed 
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that their velocity distribution was Gaussian with a dis- 
persion (including measurement error) of < 12.5 kms^^. 
Our final Draco dataset contains 112 velocities within 
39 kms~^ of Draco's mean velocity with a median ve- 
locity error of 2.4 kms~^. The union of these data 
with the dataset of IKlevna et al.l ((2001) contains 207 
unique members with good velocities. The median ve- 



locity of the combined dataset is — 290.7_ 



0.6 



kms 



where the uncertainties are obtained through bootstrap 
re-sampling. Our final UMi dataset has 144 stars within 
36 kms~^ of the mean velocity, with a median veloc- 
ity er ror of 2.9 kms~^ . Adding the earlier UMi velocities 
of Kl evna et al.l l|2003l) produces a dataset with 162 mem- 
ber stars. The median velocity of the combined dataset 
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-245.2;j;° kms" 



For both Draco and UMi, at large radii the individ- 
ual velocities of the stars relative to the mean tend to 
decline with increasing distance from the center. In the 
case of UMi, this is particularly striking with six of the 
seven outermost stars lying within la of the mean bulk 
velocity of the dSph. Figure ^ shows the radial vari- 
ation of the the line of sight velocity dispersion crp in 
Draco and UMi. The projected dispersion drops sharply 
at large radii. In each dSph, we detect no rotation be- 
yond that identified in previous work. Since our selection 
criterion for dSph membership is based on the assump- 
tion of a Gaussian velocity distribution, we would ex- 
pect to discard less than one genuine member in Draco 
(UMi) by imposing a 39 kms^^ (36 kms""'^) cut-off with 
the possible exception of extreme binaries. The presence 
of some outliers just inside the velocity cut-off at large 
radii suggests that we may, in fact, be retaining some 
non-members in our sample. If 3 percent of our data be- 
long to the Galactic foreground (assumed to have a flat 
velocity distribution between our velocity limits), then 
the overall dispersion of our sample would be increased 
by about 1 kms""'^. Thus, non-member contamination at 
large radii would tend to strengthen our conclusion of a 
falling velocity dispersion. It could be invalidated only if 
wc had erroneously removed large velocity, bound stars 
from the sample. 
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Fig. 1. — Line of sight velocity dispersion profiles (with 1-cr error 
bars) for Draco and Ursa Minor. See text for a detailed discussion. 



2.2. Photometry 

Figure 13 shows the azimuthaUy averaged surface 
brightness profile of Draco based on deep imag- 
ing to V«25 and i' «24 with the Isaac New- 
ton Telescope. This has been corrected for the ef- 
fec ts of variable extinction usin g the reddening map 
of iSchlegel. Fi nkbeiner. fc Davis! l)1998() . The hmit- 
ing magnitude of our observations is about two mag- 
nitudes fainter than that of the Sloan Digi tal Sky 
Survey data used by lOdenkirchen et aLl (|200lD to de- 
termine the light profile of Draco and permits more 
robust background subtraction in the o uter regions. 
In contrast to the Odcnkirchcn et al.l photometry, 
the profile in Figure [21 shows a clear break at ~ 
25'. We note that a similar result has re cently 
been claimed bv iKuhn. Kocevski. fc Fleckl ^2004) . The 
light profile of UMi displays a similar feature (e.g. 
Erwin & Hatzidimitriou 1995; Palma et al. 200^ at 
about 34'. This militates against the idea t hat the dSphs 
currently possess extended halos (e.g. IStoehr et al.l 
1102). 

Using simulations, IJohnston et al.l l|1999l hereinafter 
J99) suggested that stellar systems in the Milky Way 
halo will show breaks in their surface density profiles, 
beyond which unbound or extra-tidal stars begin to pre- 
dominate over bound stars. Their results are super- 
ficially similar to our photometric data on Draco and 
UMi, but there are some important differences. First, 
J99 found an enhanced velocity dispersion due to the 
extra-tidal stars and second they argued for a slow fall- 
off (oc R~^) in surface brightness bey ond the break. 
IJohnston. Choi, fc Guhathakurtal l)2002f) considered non- 
circular dSph orbits and found a wide variety of profile 
shapes and outer fall-off rates, with milder breaks and 
steeper fall-offs occurring around the apocenters of more 
eccentric orbits. They also found that the ratio of the 
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Fig. 2. — Azimuthally-averaged surface brightness profile of 
Draco. The solid and broken curves show, respectively, the best- 
fitting Plummer profile and a IKind 1196211 profile fit to the data 
within 25'. 



break radius to the actual tidal radius varies with or- 
bital phase and eccentricity and is significantly below 
unity near apocenter. ^^^ ^^^ 

Also shown in Figure|21are the best fitting lKind l|1962fl 
and Plummer profiles for the surface brightness profile of 
Draco. The former fits the inner parts well, but requires 
an additional extra-tidal population at i? > 25' to mimic 
the break. The latter provides a reasonable description 
of the entire profile. In the rest of the Letter, we use 
these models to try to understand the surprising data on 
the velocity dispersions of Draco and UMi. 



3. MODELLING 



3.1. 



Equilibrium Models 

Using the simplifying assumptions of virial equilibrium 
and spherical symmetry, the observable line of sight ve- 
locity dispersion crp as a function of projected radius R 
is 
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where I{R) is the surface brightness and v{r) is the stel- 
lar luminosity density. This expression also involves the 
mass profile M(r) of the dark matter halo, and the stellar 
velocity anisotropy parameter (3{r) = 1 — {vg)/{v^). The 
function f{r) is the integrating factor for the spherical 
Jeans equation, namely exp[— L dr2f3{r)/r]. 

Under the assumption of isotropy (/3 = 0), eqn |^ 
becomes an Abel integral equation, which can be in- 
vert ed_to_give_t he mass profile M{r) of the dark matter 
halo (iBinnev fc Trcmainc 1987, §4.2). Using an analytic 
fit to Draco's projected dispersion together with either 
a Plummer or a King profile in the Abel inversion, we 
find that the cumulative mass M{r) becomes unphysical 
{dM /dr < 0) beyond r ^ 30' . An isotropic model with a 
Plummer or King profile cannot reproduce the observed 
sharp decline in crp for Draco. Analogous fits to the dis- 
persion profile and luminosity density of UMi lead to a 
similar conclusion for its mass profile at the radius where 
its velocity dispersion falls. 

3.1.1. A Sharp Change in the Velocity Anisotropy? 
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approaching that edge. Let us suppose that as r — > rt, 
the stellar density behaves like p ^ {rt — r)". Then, by 
expanding eqn (Q), it can be shown that 
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Fig. 3. — Projected velocity dispersion for a Plummer light profile 
(scale length Rq) in a dark matter halo assuming 4 different velocity 
anisotropy laws. Even an extremely sharp change in the velocity 
anisotropy does not reproduce the decrease in anisotropy seen in 
UMi. 



One possibility is that the velocity anisotropy changes 
abruptly from isotropy (/? — 0) in the inner parts to 
strong radial anisotropy (/3 — ^ 1) in the outer parts. 
This could cause a sharp drop in the projected disper- 
sion, even if the stellar density and dark matter profile 
vary slowly and smoothly. We have investigated this op- 
tion using eqn Q together with an anisotropy parameter 
f3{r) = r"/{r"' +r2), which tends to a Heaviside function 
H{rit) about the anisotropy radius ra as n ^ cx). Figure^ 
shows the projected velocity dispersion crp assuming a 
Plummer model for the dSph luminosity density. The 
ha lo is a 'parabo lic velocity curve' model as advocated 
bv lStoehn l)2004|) as an excellent fit to sub-halos in high 
resolution simulations of the Milky Way halo. The ob- 
servable dispersion is computed under the assumptions 
of isotropy, extreme tangential or radial anisotropy to- 
gether with the Heaviside step function anisotropy. A 
very sharp change in anisotropy can cause crp to fall at 
R r^ Ti^. In fact, for a stellar population with luminos- 
ity density v oc r^" in a dark halo with an underlying 
rotation curve behaving like V{r) ~ r"''', crp falls by a 
factor of \l \/n + 27 — 2 after a sudden change to radial 
anisotropy. 

For Draco, a smooth Plummer profile provides a rea- 
sonable fit to the stellar density. Given the large error bar 
on the outermost data point of the velocity dispersion, 
a sudden change in the anisotropy at large radii could 
plausibly explain the data. For UMi, however, the out- 
ermost dispersion data point is better constrained and 
the photometry is less well represented by a smooth pro- 
file. A sharp change in the anisotropy cannot by itself 
produce the drop in crp, unless the stellar density or the 
dark halo also changes at the anisotropy radius r^. 

3.1.2. A Sharp Edge in the Light Distribution? 

If a dSph has a sharp edge in its light distribution, 
then Up would necessarily fall to zero at projected radii 



ap ^ C(rt - i?)i/2i?-(i^)/2. 



(2) 



The anisotropy (3 and the fall-off in the stellar density 
n alter the constant C, but not the scaling with dis- 
tance from the edge rt . The dark halo has been assumed 
to possess an underlying rotation curve behaving like 



V{r) 



so that the case 7 = 1/2 (Keplerian ro- 



tation) corresponds to truncation of the dark halo. So, 
irrespective of whether the dark halo is extended or trun- 
cated, the velocity dispersion of the dSph must always go 
to zero like {R — rt)^/^, if the stellar distribution has a 
sharp edge. 

For both Draco and UMi, a sharp edge to the stel- 
lar distribution seems at first sight inconsistent with the 
photometry. However, as the over-plotted King profile in 
Figure [3 shows, we can associate a tidally limited model 
with Draco provided that the excess of stars at i? > 25' 
is interpreted as an extra-tidal (and possibly unbound) 
population. If so, it is surprising that the extra-tidal 
stars are kinematically colder than those in the main 
body of Draco. This situation admits two possible ex- 
planations. First, the tails could be kinematically cold, 
either intrinsically or due to projection effects. This is 
not out of the question, as discussed later in S 13.21 Here, 
we merely note that Liouville's theorem tells us that 
the phase space density is conserved in the absence of 
mixing, so that a stretching of material to form a tidal 
tail must be accompanied by a corresponding contraction 
in velocity space. Second, given the small numbers in- 
volved, is it possible that we have missed stars associated 
with the extra-tidal enhancement? This seems unlikely, 
as the density in the tails is higher than that of Draco 
at the radii sampled by our outer velocity bin. How- 
ever, the presence of age or metallicity gradients could 
potentially lead us to sample preferentially one popula- 
tion in the outer parts - our Draco photometric data 
provide some evidence that the blue and red horizon- 
tal branches have different spatial distributions (see also 
iKlessen. Grebel. fc Harbeckll2003(l . In the case of UMi, 
the narrowness of the red giant branch makes it more 
difficult to miss a dominant tail populatio n in the outer 
parts- the weak age gradients identified bv lCarrera et alJ 
l|2002|) do not manifest themselves in the RGB. One 
might also worry that the geometry of our WYFFOS 
pointings in UMi might miss the areas dominated by the 
extra-tidal population. In fact, our outer pointings were 
aligned with the major axis of the stellar distribution and 
completely sample the elliptical region with a semi-major 
axis of 50 arcmin. 

3.1.3. Two-population Models for dSphs? 

It is worth considering whether Draco and UMi might 
contain two kinematic populations: a hot, inner "bulge- 
like" component surrounded by either a "disk" or "halo" 
component. One objection to this model is that pop- 
ulations with different spatial distributions must have 
formed under different conditions and would therefore 
not be expected to display such similar stellar popula- 
tions. However, it is possible that the weak age gradi- 
ents seen in UMi and the differences in the spatial dis- 
tributions of Blue and Red Horizontal Branch stars in 



Draco indicate the presence of more than one old popu- 
lation (see, e.g.. lHarbeck et ani20niD . 

The low projected velocity dispersion in the outer parts 
of both Draco and UMi would appear to be consistent 
with the vertical velocity dispersion of a thin disk. How- 
ever, in each case the absence of an observed rotation 
signal across the face of the dSph implies that we must 
be observing the disk approximately face-on. For the 
line-of-sight component of rotation to be smaller than 
our velocity errors, any disk would have to be aligned to 
within 6° of face-on. This is in contradiction with the 
observed flattening of UMi, which suggests that any disk 
must be close to edge-on. A pressure-supported "halo" 
population might be a more plausible candidate for the 
extra-tidal population. In Draco, the isopleths of the 
stellar density distribution suggest that the extra-tidal 
population within about 50 arcmin has the same sym- 
metry as the main body of the dSph. However, the ob- 
served kinematics require that both populations be trun- 
cated in order to give projected dispersions for each pop- 
ulation that fall rapidly in the outer parts. For UMi, 
this model suffers from the additional difficulty of the 
significant flattening required for both populations. If 
this is produced by anisotropy in the velocity distribu- 
tion (rather than by the effects of external tides) it is at 
the extreme end of the range normally seen for flattened 
systems without rotation. 

3.2. Tidal Sculpting of Draco and UMi? 

A natural explanation for a break in the light distri- 
bution of a dSph is to invoke the external tidal field pro- 
duced by the Milky Way. This perturbs the outer regions 
of all Galactic satellites leading to the escape of stars 
and the formation of tidal tails (e.g. J99). However, in 
the case of Draco and UMi there are two problems as- 
sociated with this simple scenario. First, for Draco our 
mass estimate interior to the break radius in the light is 
1 X IO^Mq while for UMi it is 2 x IO^Mq. Flattening of 
the gravitational potential and velocity anisotropy causes 
uncertainty in these masses by a factor of ~ 2. Both 
dSphs seem sufficiently massive that, given any reason- 
able, assumed Milky Way halo profile, their current tidal 
radii lie outside the observed light distribution. This sug- 
gests that neither Draco nor UMi is currently experienc- 
ing tidal disturbance of its stellar distribution. Second, 
tidal effects typically result in the heating of stellar popu- 
lations and therefore the velocity dispersion of the dSphs 
is expected to rise in the region which has been influenced 
by tides, in stark contrast to the observed data. 

The difficulties might be resolvable if the Galactic or- 
bits of Draco and UMi are significantly elongated with 
peri-Galactica smaller than 20kpc. For Draco, such 



an orbit is not inconsistent with the observed space- 
motion due to the large uncertainties on the measured 
proper motion. The proper motion of UMi is currently 
bette r constrained (Schweitzer. Cudworth. & Maiewskj 
I1997J1 and appears to rule out a deeply plunging orbit. 
However, a preliminary measurement of the proper mo- 
tion based on HST data with a four-year baseline is con- 
sistent with an elongated orbit (Piatek, priv. comm.). 
During pcrigalactic passages, the tidal radii of the dSphs 
are smaller than at their present locations. If a dSph 
passes close to the disk of the Milky Way, then the rapidly 
time-varying gravitational field disturbs its outer parts. 
This generates a heated population of stars at the edge 
of the dSph with an inflated velocity dispersion. If the 
orbit is further constrained to have a perigalacticon pas- 
sage after the mild disk shock then the heated, extra-tidal 
population can drift away from the dSph. Subsequently, 
as the dSph moves out to its current location, the tidal 
radius increases again and the colder stars in the tidal 
tails which have not had sufficient time to escape during 
the pericenter passage are recaptured. This results in a 
population which has the morphological appearance of 
a tidal tail but is kinematically cold. It is also possible 
that the cold clump near the center of UMi IjKlevna et alJ 
I2003J) is a projection of this cold, extra-tidal population 
onto the face of the dSph. 

This scenario places rather tight constraints on the 
properties of both Draco and UMi as well as on their pos- 
sible Galactocentric orbits. First, their orbits must take 
them within 20 kpc of the Milky Way. Second, a simple 
estim ate of their tidal radii (following, e.g., Klcvna et aj 
|2001[) shows that even at this radius their tidal limits are 
large with respect to their stellar populations unless their 
masses are below 5 x IO^'Mq, which is at the lower limit 
allowed by the data. Third, their orbits must spend suf- 
ficient time near pericenter to allow the stars which have 
been most heated during the close passage to escape be- 
fore the sphere of influence of the dSph engulfs them. We 
conclude that this scenario is only plausible if Draco and 
UMi have somewhat lower masses than previously esti- 
mated and have a low dark matter density outside about 
30 arcmin. Further both dSphs must be on deeply plung- 
ing orbits which take them close to the disk of the Milky 
Way. We are currently performing iV-body simulations 
to investigate this scenario in more detail. 
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